Entanglement in spin-1/2 dimerized Heisenberg systems 
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We study entanglement in dimerized Heisenberg systems. In particular, we give exact results of 
ground-state pairwise entanglement for the four-qubit model by identifying a Zi symmetry. Al- 
though the entanglements cannot identify the critical point of the system, the mean entanglement 
of nearest-neighbor qubits really does, namely, it reaches a maximum at the critical point. 

PACS numbers: 03.65.Ud, 03.67.-a 
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The study of entanglement properties of many-body 
condensed matter systems has received much atten- 
tion Q-jslJ. It covers both the rapid developing area 
of quantum information and the area of condensed mat- 
ter. Various quantum models have been considered in the 
literature and some of them are exactly solvable. Specif- 
ically, it was well-known that the anisotropic Heisen- 
berg model can be solved formally by Bethe's Ansatz 
method [3ll l3^| for arbitrary number of qubits N, how- 
ever, we have to solve a set of transcendental equations. 

It is interesting to see that entanglement in a system 
with a few sites displays general features of entanglement 
with more sites. In the anisotropic Heisenberg model 
with a large number of qubits, the pairwise entanglement 
shows a maximum in the isotropic point |33| . This fea- 
ture was already shown in a small system with four or 
five qubits |34j. So, the study of small systems is mean- 
ingful in the study of entanglement as they can reflect 
general features of larger or macroscopic systems. 

Dimerized systems play an important role in condensed 
matter theory, and entanglement in dimerized systems 
has been considered [2^,|33- Here, we aim at obtaining 
some analytical results of entanglement in a small system. 
Chen et al [3i| found that the concurrence [3(|, quanti- 
fying entanglement of a pair of qubits cannot enable us 
to identify the critical point J2/J1 = 1- However, we 
will show that although the single concurrence cannot be 
used to identify the critical point, the mean concurrence 
really does. 

For N < 7, the isotropic Heisenberg Hamiltonian can 
be analytically solved [371 |3j| due to the translational 
invariant symmetry which can be used to reduce the 
Hamiltonian matrix to smaller submatrices by a factor of 
N (3^. However, for our dimerized systems, the trans- 
lational invariant symmetry breaks. We identify a Z2 
symmetry in the dimerized model, and using this sym- 
metry we can further reduce Hamiltonian matrices by a 
factor of two. We will first consider the simplest repre- 
sentative four-qubit dimerized system, and by using the 
Z 2 symmetry, all eigenvalues can be obtained. 

Eigenvalue problem. The Heisenberg Hamiltonian for 



the dimerized chain with even number of qubits N reads 



N/2 



^l(2S2i_l • S2i + — ) + J2(2S2j • S2i+1 + — ) 



H =J2 

i=l 

N/2 

= / ^(Jl«^2t-l,2j + J2S 2 i,2i+l) 



(1) 



i = l 



where Si is the spin- half operator for qubit i , <Sjj+i = 
^ (1 + <jj • <Xj+i) is the swap operator between qubit i and 
j, <7i = (<Ji X , o~iy, Oiz) is the vector of Pauli matrices, and 
Ji and J2 are the exchange constants. We have assumed 
the periodic boundary condition, i.e., TV + 1 = 1. In the 
following discussions, we also assume J\ , J2 > (antifer- 
romagnetic case). 

For the four-qubit case, the dimerized Hamiltonian 
simplifies to 



H = Ji(iSi2 + 534) + J 2 {S 2 3 + S41). 



(2) 



Although we have imposed the periodic boundary condi- 
tion, the Hamiltonian is no longer translational invariant 
except the case of J\ = Ji- As we stated above, the 
translational invariant symmetry can be used to reduce 
the Hamiltonian matrix to smaller submatrices by a fac- 
tor of N j^. However, for our case, we cannot have this 
reduction. 

To exactly solve the eigenvalue problem of the Heisen- 
berg model, we first note that [H, J z ] = 0, the whole 16- 
dimensional Hilbert space can be divided into invariant 
subspaces spanned by vectors with a fixed number of re- 
versed spins. Then, the largest subspace is 6-dimensional 
with 2 reversed spins. Here, J z — J2i=i <J iz/^- Due to 
the Z2 symmetry [H, S x ] = 0, it is sufficient to solve 
the eigenvalue problems in the subspaces with r reversed 
spins, where r E {0, 1, 2} and 



S T <Jr 



(3) 



The subspace with r = only contains one vector 
1 0000), which is the eigenvector with eigenvalue 



E = 2(Ji + J 2 



(4) 



For subspace with 2 > r > 0, the smallest submatrix 
is 4 x 4. We need further reduce submatrices to smaller 



2 



ones. Although we do not have a translational invariance, 
from the four-qubit Hamiltonian J2J , we can identify an- 
other symmetry Zi symmetry, given by 



[H, £ 



12,34J 



[H, S13S2 



0. 



(5) 



The operator exchange the state of qubits 1 and 2 with 
the state of qubits 3 and 4, namely, 

Si2,34|™-i, r7 i2,TO 3 ,m 4 ) = |m 3 ,m4,mi,m 2 ). (6) 

The eigenvalue of the operator £12,34 on all energy eigen- 
states is either 1 or -1. 

The subspace with r = 1 is spanned by four basis vec- 
tors { 1 1000), |0100), |0010), |0001)}. Taking into account 
the Z2 symmetry, we choose the following basis 



1 



i± 



-^=(11000) ± |0010», 



>i±) = -^=(10100) ± 10001)). 



(7) 



Obviously, they are all eigenstates of operator £12,34, and 
basis {|^+), |0+)} ({IV 7 -), 1^-)}) spans an invariant two- 
dimensional subspace of Hamiltonian H . Let the Hamil- 
tonian act on the basis, then we obtain 

H\tf> 1+ ) =(Ji + J 2 )(|Vi+)+0i+)), 
H\<h+) =(Jl + J2)(|Vl+)+01+)), 

#|^i_) =(Ji + J2M1-) + (Ji - J2M1-), 
tf|0i_) =(Ji + J2M1-) + (Ji ~ J 2 )|^i-). (8) 

From the above equation, the eigenvalues of Hamiltonian 
H is obtained as 

£1,1 = 2(J X + J 2 ), E h2 = 0, £7 lj3 = 2Ji, E 1<4 = 2J 2 , (9) 

where the first subscript denotes the number of reversed 
spins. 

For the case of r = 2, due to the existence of the 
Z2 symmetry, we choose the following basis for the six- 
dimensional subspace 

|V*±> = ^(liioo) ± loon)), 



1 
1 

71 



(|1001) ± |0110)), 
(|0101) ± |1010)). 



(10) 



Then, we can reduce the 6x6 Hamiltonian matrix to a 
block-diagonal form with two 3x3 matrices. After the 
action of the Hamiltonian on the above basis, we obtain 

H\i; 2 _) =2Ji|V> 2 _), 
H|0a_) =2J 2 \h-), 
H\<P2-) -0, 

H\i)2+) =2Ji|^ 2 +) +2J 2 |(^2+), 
ff|02+) =2J 2 |02+) +2Ji|^ 2+ ), 

H\<p 2+ ) =2J 2 |V> 2 +) + 2Ji|0 2+ ). (11) 



One 3x3 block is of the diagonal form, and the eigen- 
values read 

Ea,i = 2Ji, E2.2 = 2J 2 , £2,3 = 0. (12) 

Another block is written as 

2Ji 2J 2 
H= I 2J 2 2 Jx ) . (13) 
2J 2 2Ji 

Although this is a 3 x 3 matrix, we can further reduce 
it to 2 x 2 matrix since, from the last three equations 
of Eq. (|llfl . one eigenvalue is found to be £2,4 = 2(Ji + 
J2). With the help is this eigenvalue, the characteristic 
polynomial of the Hamiltonian matrix can be brought to 
a quadratic form by dividing it with E — -£?2,4- Then, we 
obtain 



£2,5 =2^ J( + Ji - J1J2 

EgS =E 2 r 



2JJ( + Jj- JxJ 2 . 



(14) 



Thus, all eigenvalues are analytically obtained for the 
four-spin dimerized Heisenberg model. We see that 
ground state is non-degenerate and the energy is given 
by £2,6- Although the eigenstates can be easily obtained, 
they are not given explicitly here as the knowledge of 
eigenvalues is sufficient for discussions of entanglement 
properties as we will see below. 

Entanglement in the four-qubit model. We first study 
the ground-state entanglement of qubits 1 and 2. Due 
to the SU(2) symmetry in our Hamiltonian, the concur- 
rence quantifying the entanglement of two qubits is given 

by m 



a 



-2(Si ■ S 



i+l) 



l/2 = -(£ 



M+i/' 



(15) 



where we have ignored the max function in the usual 
definition of the concurrence, and thus the negative con- 
currence implies no entanglement. We see that the en- 
tanglement is determined by the expectation value of the 
SWAP operator. There are two distinctive concurrences 
C12 and C23 , corresponding to two spins coupled by bond 
J 1 and J2. 

The expectation value {S12) can be calculated via 
Feynman-Hellman theorem. By applying the theorem 
to the ground state, we obtain 



{S12) + (£34) 



dE, 



GS 



J 2 - 2 Jx 



9Ji VjJ 



J1J2 



(16) 



Due to the Z2 symmetry [H, £12,34], we have (£12) = 
(£34). Then, from the above equation, the expectation 
value (£12) is obtained. Substituting it to Eq. (|T5)l leads 
to 



Cxi — 



2Jx - J 2 



2^1f 



J 2 



J1J2 



(17) 



FIG. 1: The concurrences versus J2 in the four-qubit dimer- 
ized model (Ji = 1). 



FIG. 2: The concurrences versus J2 in the six-qubit dimerized 
model (Ji = 1). 



In a similar way, we can obtain the concurrence C23 = 
C14 as 



2J 2 - Jj 



2:( 



(18) 



Thus, the analytical expressions of two types of concur- 
rence are obtained, and the concurrence Ci 2 can be trans- 
formed to C23 by exchanging J\ and J 2 . 
From Eqs. ljT7|) and (JTSJl, we read that 



= 1 when J 2 = 0, 
C12 4 = !/ 2 when J 2 = Ji, 
< when J 2 > J, 



(12) 
2th 



(19) 



2Ji, 



and 



< when 3 2 < jg® = J x /2, 
C23 { — 1/2 when J 2 = Ji, 
= 1 when J2 = 00, 



(20) 



In the case of J2 = 0, the ground state is just the un- 
coupled two singlets, and thus C\i = 1. When J 2 = Ji, 
the dimerized Hamiltonian is reduced to the isotropic 

one, and the concurrence C12 = C23 = 1/2 6]. When 

(23) 

J 2 > Jjjth > ^ nc entanglement between qubits 2 and 3 

(12) 

builds up, and when J2 > J 2th = , the entanglement 
between qubits 1 and 2 vanishes. 

It was argued that the pairwise entanglements quanti- 
fied by the concurrence cannot ide ntify the critical point 
(J 2 = Ji) of the dimerized system |35| . however here we 
use the mean entanglement of nearest-neighbor qubits as 
the system we are considering is not a uniform system, 
namely, it has two distinct exchange constants. There 
are four pairs of nearest-neighbor qubits, and thus 

^ _ C12 + C23 + C34 + C41 _ C12 + C 2 3 . . 



Substituting Eqs. lfT7|l and (fT5f> to the above equation 
leads to 



from which we read 



J 1 + J 2 



(22) 



1/4 when J 2 = 0, 
Cmoan = { 1/2 when J 2 = Ji, 
1/4 when J2 = 00. 



(23) 



In Fig. 1, we plot the concurrences and the mean con- 
currence versus J2- We see that the mean entanglement 
takes its maximum at the critical point J 2 =Ji. From 
Eq. (|22|l . we can also check that the maximal point is at 
the critical point. In the following, we will numerically 
show that this fact still holds for dimerized chains with 
more qubits. 

Numerical results. Now we consider more general cases 
of even-number qubits. For even N > 4 , the analytical 
results of entanglement are hard to obtain, and here, we 
use exact diagonalization method to numerically calcu- 
late the entanglement. For even-number case, the mean 
entanglement is still given by Eq. (|21|) . 

In Figs. 2 and 3, we plot the concurrences versus J2 
for N — 6 and N = 8, respectively We observe that the 
entanglement properties are similar to those in the four- 
qubit model, namely, they are qualitatively the same. 
The mean entanglement reaches its maximum at the crit- 
ical point. From Figs. 1-3, we see that the threshold value 



^2th decreases, while J 2 ^ ; increases when the number of 
qubits increases. 

Finally, we give a relation between the two types of 
entanglement with the ground-state energy Eqs. From 
Eq. we immediately have the relation between the 
ground-state energy and the correlators (<Si 2 ) and (5 2 3), 



,-(23) • 

/ V ' IT 



4 




FIG. 3: The concurrences versus Ji in the eight-qubit dimer- 
ized model (Ji = 1). 



E G s/N = -(Ji(5 12 ) + J 2 (S 23 ». (24) 
Then, from Eq. I|15|) . we obtain 

Egs/N = -\{JiC 12 + J 2 C 23 ). (25) 
The above equation gives the relation between ground- 



state energy and the two typical concurrences. If we know 
the entanglement of qubits 1 and 2 and the ground-state 
energy, we can know the entanglement of qubits 2 and 3 
from the relation. 

Conclusion. We have studied entanglement in spin- 
1/2 dimerized Heisenberg systems. As a representative 
system, the four-qubit model was studied in detail. By 
identifying the Z 2 symmetry, we have solved the eigen- 
value problem completely. From the ground-state energy, 
the analytical results of the two types of pairwise entan- 
glement were obtained. We have located the threshold 
value of J 2 after which the entanglement between qubits 
1 and 2 vanishes and the threshold value before which 
there exists no entanglement between qubits 2 and 3. 

In order to identify the critical point of the dimer- 
ized system, we have proposed the mean entanglement 
of nearest-neighbor qubits as an efficient indicator. The 
mean entanglement displays a maximum at the critical 
point. We have also considered the case of more qubits, 
and the numerical results show that the entanglement 
properties are similar to those in the four-qubit model. 
Dimerized Heisenberg systems play an important role in 
condensed matter physics. It is interesting to study en- 
tanglement in other dimerized systems, which is under 
consideration. 
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